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Abstract 

We obtain an algorithm computing explicitly the values of the non solvable spectral radii 
of convergence of the solutions of a differential module over a point of type 2, 3 or 4 of 
the Berkovich affine line. 

Introduction 

By a theorem of Young [You92] the small values of the radii of convergence of the solutions of a 
differential operator are explicit, and coincides with the small slopes of the Newton polygon of a 
differential operator attached to the module (cf. Prop. 2. 2). Larger radii are not immediately readable 
on the coefficients of the operator. This discrepancy is a peculiarity of the p-adic world: only small 
radii are " visible" . To overcome this problem B.Dwork observed that the pull-back by Frobenius 
functor increase the radii of the solutions, and together with G.Christol [CD94] they constructed 
an inverse of the Frobenius functor (often called Frobenius antecedent) in order to make the radii 
of the solutions smaller, and hence explicitly intelligible in a cyclic basis. Although theoretically 
satisfactory, the inversion of Frobenius is a completely implicit operation. Moreover the antecedent 
exists only if all the radii of the solutions are not small. So one is obliged to factorize the module 
by the radii of the solutions if one wants to understand the non minimal radii of the solutions. 
The factorization is also an implicit operation. Recently in [KedlO] K.Kedlaya observed that the 
Frobenius Push-forward operation has essentially the same effects as the inversion of the Frobenius 
on the radii of the solutions, and he is able to control the exact behavior of all the radii of the 
solutions under this operation (even small radii). 1 Frobenius push-forward functor is completely 
explicit, and it allows to obtain a concrete algorithm to compute the radii of the solutions that 
are not maximal (i.e. non solvable). The price to pay is that the dimension of the push-forward by 
Frobenius is p-times that of the original module, so that the complexity of the algorithm is multiplied 
by p at each application of the push-forward. 2 Hence more the radius is large, more the complexity 
increase. Moreover the algorithm admits an end if and only if all the radii are not maximal (i.e. non 
solvable). Eventually we provide the algorithm, but we avoid to provide a complete formula as one 
does for rank one equations (cf. [Chrll] 3 ). Indeed the complexity seems so great that the formula 
would result not useful to be written. Explicit examples are quite complicate even in the rank one 
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The very first reference for this is [Chr77], in which the author introduces the push-forward, and its relation with 
the pull-back, and use it to prove the inversion of Frobenius (existence of the antecedent) under the condition that 
the solutions of the module are Bounded. 

2 Xavier Caruso recently pointed out that the explicit factorization of an operator by the radii of convergence seems to 
be concretely implementable into a machine. This would highly reduce the complexity of the present algorithm since 
by considering the right factor of the push-forward by Frobenius the dimension remains constant at each step. 
3 The formula that we have contributed to prove in [Chrll] is based on a completely different approach, and it uses 
Witt vectors (following techniques of [Pul07]) to explicitly describe Taylor solutions of a rank one differential equation. 
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case. This note is intended to make explicit the computations, and the link between the different 
results, in view to make them explicitly calculable by a computer. 

0.0.1 Acknowledgements. We thanks Kiran S. Kedlaya for Remark 3.3. 

1. Radii of convergence 

Let (K, |.|) be a complete field with respect to an ultrametric absolute value |.| : K — > R^o- Let L/K 
be a complete valued field extension, let c G L, and p > 0. For all polynomial P(T) := ^ cijT* G 
K[T] define |P| CjP := sup n>0 \P® (c)\p n /nl The setting \Pi/P 2 \ c , P := \Pi\c, P I \P2\cp defines an ab- 
solute value on the field of fraction K(T) := Fiac(K[T]), and hence a Berkovich point £ C)/9 of the 
affine line A^ an . Since p > if c G K one obtains in this way all the points of type 2 or 3 of A^- an , if 
one allows c ^ K one also has all points of type 4. The derivative d/dT extends by continuity to the 
completion J$f c ,p of (K(T), |-| C p)- A differential module over J4? c ,p is a finite dimensional J4? c ^ p - vector 
space M together with a iT-linear map V : M — > M satisfying V(/m) = d(f)m + /V(m), for all 
/ G J^, p , m G M. Let r(£ C)P ) ^ p be the radius of the point £ CjP (cf. [Pull2, Section 1.3.1]). If Q/K 
is a complete valued field extension such that Ao' an has an il-rational point i CjP G SI lifting £ CjP , then 
r (£c,p) is the radius of the largest open disk DQ(t CjP , r(£ CiP )) satisfying D^(i cp , r(£ CjP )) n i^ alg = 0. 

Definition 1.1. Let r := rank(M). For i = l,...,r we denote by Hi = 7£f' sp (£ c , P ) ^ r(£ CjP ) 
the radius of the largest open disk in D^(t CjP , r(£ CjP )) centered at t CjP on which M has at ieast 
r — i + 1 Vt-linearly independent Taylor solutions (cf. [Pull2, section 4.2] or [KedlO, 11.9]). One has 
Tlx ^ H 2 < • • • < Ur- 

We say that Hi is solvable if TZi = r(£ cp ). In this paper we provide an algorithm computing non 
solvable radii. 

2. Comparison of Newton polygons and computation of small radii 

Let r ^ 1 be an integer. A slope sequence is the data of r real numbers s\ ^ . . . ^ s r in increasing 
order. Define the i-th partial height as hi := si + • • • + Sj. A slope sequence defines univocally a 
convex function h : [0,r] — ¥ M by /i(0) := 0, h(i) := /ij, and h(x) = SiX + (hi — i-Si) for all x — l,i], 
i = 0, . . . , r. The function /i is called the Newton polygon with slopes s% ^ . . . ^ s r . 

Definition 2.1. The Newton polygon with slopes Sj := sf ,sp (£ C)P ) := m(7Zf ,sp (£ c ,p)) is cabled the 
spectral Newton polygon of M. We denote by hi := /if' sp (£ c ,p) its i-th partial height. 

Let C = Yli=o 9r-id l , G ^ iP , be a differential operator with go = 1 and fly 7^ 0. Let vo = 0, 
and for all i = l,...,r let Vi := — ln(|<7j| CiP /cj l ), where w := lim n | n! 1 1 / ri . Let Li := {(x,y) G 
M? such that x = i,y ^ Vi}, note that Lj is empty if and only if gi = 0. Define the spectral Newton 
polygon NP(C) as the intersection of all upper half planes H a ^ := {(x,y) G M 2 such that y ^ 
ax + 6} with Vj)}i=o,...,r C -ff a ,fe- Let /i £ : [0, r] — )■ E be the convex function whose epigraph is 
NP(C): h c (x) = min{y such that (x, y) G NP(C)}. Explicitly one has hf := h (i) = sup sgK {s • i + 
minj = Q^ ^ r Vj — s ■ j}. Then NP(C) is the Newton polygon with slopes {sf := hf — hf_ 1 }i=i r .. jr . 

Proposition 2.2 ([You92]). Let £ be a differential operator as above and let M be the differential 
module dehned by C. Let C := ln(cj • r(£ CjP )), then for all i = 1, . . . , r one has 

min(sf ' sp , C) = min(sf,C) . □ (2.1) 

Remark 2.3. In order to apply (2.1) we need an algorithm to find a cyclic basis of M (cf. section 
3). If the absolute value of K is trivial on Z (i.e. if \n\ = 1 for all n G Z — {0}), then uj = 1, and 
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Proposition 2.2 allows to find all the radii 72-j. If the absolute value of K is p-adic (i.e. if \p\ < 1), 



then uj = < 1. In this case, we also need a technique (Frohenius push-forward) making the 

(non solvable) radii smaller than u:r(^ CyP ) (cf. section 4). 



Let (F,d) be a differential field and let F(d) = (Bi^oF o d % be the Weil algebra of differential 
operators. The multiplication of F(d) extends that of F = F o d° by the rule do f = f o d + d(f), 
for all / 6 F. Finite dimensional differential modules over F are exactly torsion left F(d) -modules. 
The so called cyclic vector theorem asserts that all differential module are not only torsion modules 
over F(d) : they are cyclic modules i.e. of the form (M, V) = (F{d)/F{d)C,d), for some C := 
Y^4=o9r-id l £ F{d), with go = 1, gi £ F. The image of {l,d,d 2 , . . . ,d r ~ 1 } in the quotient form a 
basis of M, and the action of V is given by the multiplication by d in the quotient. In fact the cyclic 
vector theorem is equivalent to the existence of an element c G M, called cyclic vector, such that 
{c,V(c),V 2 (c),... ,V r-1 (c)} is a basis of M. In this case if cj := V*(c), and if V r (c ) = Ei=o 
then fi = —g r -i. The existence of such a vector is due to [Del70, Ch.II, Lemme 1.3]. Subsequently 
N.M.Katz provided the following explicit algorithm 

Theorem 3.1 ([Kat87]). Let (M, V) be a differential module over (F,d) of rank r, and let e := 
{eo, • • • , e r _i} C M be a basis of M. Let do, . . . , a r ( r -x) £ F be r(r — 1) + 1 distinct constants i.e. 
d(o>i) = 0. Then at least one of the following elements of M is a cyclic vector: 



Remark 3.2. The explicit base change matrices to the Katz's cyclic basis are quite complicate and 
have been explicitly computed in [Pul]. On the other hand the proof of the existence of a cyclic 
vector of Deligne [Del70, Ch.II, Lemme 1.3] proves that the family of cyclic vectors in a given 
module is the complement of an hypersurface. The base change matrices of Katz's algorithm are 
quite involved and hard to hnd by hand even in small examples. It is often convenient to pick an 
arbitrary vector and test if it is cyclic. 

Remark 3.3. One shall avoid the use of a cyclic basis using [KedlO, Lemma 6.7.3, Thm.6.7.4, 
Conjecture 4.4.9]. 



In this section we assume that \p\ < 1 (cf. Remark 2.3). 

Hypothesis 4.1. TzY' sp is insensitive to scalar extensions of K, and by translations. So in the 
sequel we will assume c = and replace the indexation (c,p) by p. In this case one has r(£ p ) = p. 
We then work with \.\ p , £ p , Jtf? p , r(£„) = p with the evident meaning of notation. If the reader 
needs to preserve the setting (c,p), the same computations hold replacing the map <p : T h-> T p by 
T i — y (T — c) p + c. Or alternatively one also can preserve (p : T — > T p , and proceed as in [Pull2, 
Section 7] to check the behavior of the radii by Frobenius at points that are close enough to the 
principal branch p^r \-\o p (this is often necessary if one search the slopes of 7^' sp along a Berkovich 
path p i y \-\c,p, with c € K and p close to \c\). 

Let T,T be two variables, and let (p : K(T) — > K(T) be the ring morphism sending T into T p . 
This extends into a isometric inclusion tp : J%? pP — > J^f p of degree p. One has the rule ^(/(T)) = 



3. Explicit Cyclic vector 




□ 



(3.1) 



4. Frobenius Push-Forward and explicit computation of larger radii 
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0L(J(T)), for all / G K(T). We call 
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dT 



(4.1) 



Let (M, V) is a differential module over (J^ p , d p ) of rank r. Since {d p v)\ x p = d pP , then (M, (pT p l ) X V) 



is a differential module over (JUp, d p p) that can be seen (by restriction of the scalars) as a differential 
module over (J#pp,d p p) of rank pr. We call (<£>*M, <^*V) the differential module so obtained. 



4.1 Explicit matrix of <^*V. 

One has a direct sum decomposition J4? p ■ 

uniquely written as g{T) = 5Zfc=o 9k{T p )T k 



® p k=Q ip{M'pp) ■ T k , so that each g(T) G J^f p can be 
Y^k=o 9k{T)T k . The derivation d pP stabilizes globally 



each factor and d p v{gk{T)T ) = d p p{gk{T))T . For all g(T) G Jt? p we define (p*(g)(T) € M pxp {M ? p p) 
as the matrix of the multiplication by 5(T)/(pT p_1 ), with respect to the basis 1, T, ... , T p_1 over 
JtfpP. One has 



/ g p -i{T) Tg p _ 2 (T) Tg p _ 3 (T) 
I g (T) g p -i(T) Tg p ^ 2 (T) Tg p - Z (T) 



v., 



TgiiT) v 
Tg2(T) 



(4.2) 



gi(T) g (T) g p -i(T) Tg p . 2 (T) Tg p ^{T) 
.g P -2(T) g p -z{T) - - S i(T) 90 (T) Sp-i(T)/ 

Notice that the terms over the diagonal are multiplied by T. Let (M, V) be a differential module 
over J4? p . Fix a ,J%? p -lme&r isomorphism Jt? p r 

V in this basis, where G(T) = (gij^T))i,j=i,...,r G M rxr (Jif p ). Writing Jif p 1 
one sees that the matrix of </?*(V) is given by the block matrix 

ip*{G)(T) := (</'*(5i,i)(r))i,j=i 1 ...,r G M prxpr (J^ p p) . (4.3) 



M (i.e. a basis of M), and let -i — G(T) be the map 

" {® p k 7M^ P p)-f k r 



4.2 Behavior of the radii by Frobenius push- forward 

Theorem 4.2 ([KedlO, Thm.10.5.1]). Let fti < . . . < TZ r be the radii of M at £ p (cf. Def.1.1). Let 
i\ be such that TZi x ^ up < IZ^+i. Then, up to permutation, the spectral radii of <^*M at £ p j> are 

J { h^- 1 ^,...,!^- 1 ^ } (J (^, ^v,...^v) . □ (4.4) 

p-times * >n p-l-times 

If sx ^ • • • ^ s r is the slope sequence of the spectral Newton polygon of M at £ p , and if io ^ i\ 
satisfies 1Zi < p = 7Zi 0+ i, 5 then by Theorem (4.2) the slope sequence associated to (p*M at £ p p is 

p-times p-times 



InfH/- 1 ) + Sl = ■ ■ ■ = lnfjpl/- 1 ) — < ln(\p\ff-^ + s h = ■ ■ ■ = In^f?' 1 ) + s n < 

(p— l)(r— ii)-times , ■ w 
Ke ' K ' (r— io)-times 

^ ln(wV) = • • • = ln(wV) < ps il+1 ^ • • • ^ ps io < ln(/) = • • • = ln(/) . (4.5) 

We have two main goals here. Firstly the sequence s\ ^ • • • ^ s r is perfectly determined by the 
knowledge of the slope sequence (4.5) of (p*M (even if some of the Sj are equal to the critical value 
ln(w/o)), see [Pull2, Prop. 6.17] for a more precise statement. Secondly the values of s, satisfying 

4 It is understood that ii = if ojp < IZi. 
It is understood that io = r if 1Z r < p. 
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ln(w 1 / p / o) ^ Si < ln(w 1 / p / o) corresponds to small radii 6 of tp*M that are explicitly intelligible by 
Proposition 2.2. Iterating this construction by performing several times the push-forward one obtains 
an explicit algorithm that computes all the non solvable radii TZi, . . . , 1Zi in a finite number of steps. 
Once this have been achieved, one knows in fact all the spectral radii since the remaining radii are 
all equal to p. Unfortunately Proposition 2.2 does not furnish any information about radii that are 
larger that up, so (unless the radii are all not solvable) it seems impossible to know whether the 
algorithm is ended or if one needs more applications of the Frobeinius push-forward. 
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